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Abstract 

A property of dynamical correlation functions for nonequilibrium states is discussed. We 
consider arbitrary dimensional quantum spin systems with local interaction and trans- 
lationally invariant states with nonvanishing current over them. A correlation function 
between local charge and local Hamiltonian at different spacetime points is shown to 
exhibit slow decay. 
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1 Introduction 

In this paper we treat nonequilibrium states of quantum spin systems. When we speak of 
nonequihbrium states, there are two kinds of states, that is, states near equihbrium and 
far from equihbrium. As for the investigation of the former nonequihbrium states, the 
main purpose is to understand how the states approach to an equihbrium state. The latter 
arc states which, on the other hand, can not be treated as perturbed equihbrium states 
and do not converge to any equihbrium states, whose properties arc less known. Wc, in 
this paper, study a property of the latter nonequilibrium states. There have been works on 
the nonequilibrium states from various points of view [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. Among 
the various results, it is known that nonequilibrium steady states exhibit slow decay of 
equal-time correlation function from approximate theories like fluctuating hydrodynamics 

[I] . Recently an approach to the problem from exactly solvable models [2] has been 
investigated and the expected behavior of equal-time correlation function was confirmed. 

In this paper, another type of slow decay of dynamical correlation function will be 
discussed. We consider locally interacting quantum spin systems and translationally in- 
variant states with nonvanishing current over them. We do not impose any other con- 
ditions like stationarity of the states and do not ask how the nonequilibrium states are 
realized. Since the states have finite current at infinity, they cannot be locally perturbed 
equilibrium states and thus are far from equilibrium. We consider a correlation function 
between local charge and local Hamiltonian at different spacetime points. It will be shown 
in mathematically rigorous way that the correlation functions of the states exhibit slow 
decay with respect to space and time. The way to show the claim is a generalization of 
the method employed in [10]. 

The paper is organized as follows. In section 2, we introduce our nonequilibrium states 
on quantum spin systems. In section 3, the main theorem is proved. 

2 States with Nonvanishing Current on Quantum 
Spin Systems 

Let us consider a d-dimensional {d — 1, 2, 3, • • •) infinite quantum spin system (see e.g., 

[II] ). For simplicity, we assume nearest neighbor translationally invariant interaction. 
That is, for each pair of neighboring sites x and y (i.e.,x, y G Z'^ s.t. |x — y| = 1), a 
self adjoint operator /i(x, y) is defined and satisfies Tz (/i(x, y)) = /i(x + z, y -h z) for all 
z e Z*^ where stands for the translation operator. Denoting fundamental vector with 
direction i as ej, local Hamiltonian with respect to a finite region A C Z is defined as 

d 

r=l {x,x+eT.}cA 
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In this paper we employ the Heisenberg picture and the Hamiltonian defines time evolu- 
tion of local operators, say A, by 

at(A) := lim e'^^^Ae''^^^ 

where the limit is taken with respect to norm topology (we put h = 1). To discuss 
(electric) current, a charge n(x), a self-adjoint operator, is assumed to be defined on 
each site x [12], and the charge on a finite region A is denoted as A^a '■= Z]xeA^(^)- 
It induces a symmetry transformation on each finite sublattice A, that is, [Ha, Na\ — 
holds. The current operator between sites x and x -|- is defined by j(x, x -|- e^) := 
i[n(x), /i(x, x + e,.)] = —i[n{-x + e^), /i(x, x + e^)]. 

The above seemingly abstract setting has physically interesting examples. For instance, 
interacting fermion system is on the list. For each x G Z"*, charge is defined as n(x) := c*Cx 
where Cx is a fermionic annihilation operator at site x and /i(x, x -|- e^) = — T(c*_,_e^Cx -l- 
CxCx+e^) -l- Vn{'x.)n{'x + er) gives a nearest neighbor Hamiltonian where T and V are real 
constants. The current operator is calculated as j(x, x + e^) = iT{c^_^_^^Cx — c*Cx+e^)- 
The Heisenberg model can be another example. On each site x G Z, spin operators 
(^S^\ S^''^ satisfying usual commutation relations for angular momentum live on. 
/i(x, X + e^) := 4^)S'Se, + S^^S^^l^^ + A5(3)^(3)^^ ig ^ ^g^^ constant) and n(x) := 4^) 
leads the current j(x, x + e,) = -S^^S^%^ + S^^^S^^l^^. 

In this paper, we consider translationally invariant states with nonvanishing current at 
time t = 0. Without loss of generality, wc can fix the first axis ei as a direction of 
nonvanishing current. That is, we consider translationally invariant states which satisfy 

0-(O,ei))^0. 

Since our conditions are weak, it surely contains a several physically interesting models 
[3, 4]. However, we do not impose any other conditions like stationarity and stability, 
and therefore some states might be hardly realized in real physical situation [13]. The 
crucial point is that the states have nonvanishing current at infinity and they are not, 
say, locally perturbed equilibrium states. 

3 Slow decay of correlation function 

In derivation of slow decay of correlation functions for continuous symmetry breaking 
equilibrium states (Goldstonc theorem), it is crucial to represent an order parameter by 
a commutator [14, 15]. Our discussion can be considered as a nonequilibrium analogue 
of [14]. For nonequilibrium states, the following is an important observation. For finite 
sublattices V and A satisfying V D A, 

[-i ^ a;in(x),i/A] = ^ j(y,y + ei) (1) 

XGV {y,y+ei}CA 

holds where xi is a component of x = {xi,X2, ■ ■ ■ ,Xd) G Z'^. Hereafter, all the finite 
sublattices are assumed to be d-dimensional cubes, and we use simplified notations Hr :— 



3 



H[-R,R]d for positive integer R where [—R, R^ '■= {x = {xi, X2, - ■ ■ , Xd) & 'L'''\ — R < Xj < 
R for all j — 1,2, ■ ■ ■ , d}. Taking expectation value of (1) with respect to a translationally 
invariant state and we put A = [—R, R]'^, take W > R and divide (1) by |{y|{y, y + ei} C 
[-R,R]'^}\ = 2R{2R+l)'^-\ Then we obtain 

-— ^(H Yl xin{^),HR]) = {j{0,ei)), 

where we have used the translational invariance of the state. By letting the size of cubic 
lattices in the above equation infinity, we obtain 

lim lim — ^- — -rT(H Y ^x), //i?]) = (j(0, ei)). (2) 

R^^w^^2R{2R+lY-^^^ xehW 

The ordering of the above limiting procedures is crucial and can not be exchanged. In 
fact it is easy to see that if we take i? — > oo first, for steady states it always gives zero. 
Moreover we can show the following lemma. 

Lemma 1 For an arbitrary translationally invariant state and for any time t e R, 
op/op^^iW-i ^[~^ ^ xin(x),Q;t(i/fl)]) = (j(0,ei)) 



lim lim 

R—*oo W-^oo 



holds, where at{HR) is a time evolved object of Hr (remind that we are in the Heisenberg 
picture). Although the Ihs includes time t explicitly, its value does not depend on t. 

Proof: 

Since for t = 0, the above equation holds (see (2)), we estimate its deviation for finite t. 
That is, since 

\-i xin{yi),at{HR)] 



2R(2R + ir-^^ 



[-i xin{x),HR] 



2R{2R + ir-^^ xehw" 

+ ogrop^iw-i / Yl xin{x),-ias{[HR,HR+i])] 

znyzn -h ij JO ^e[-w,w\d 

holds, to show 

lim lim op.op^, TW_i llH Y x^n{^),-ias{[HR,HR+^\)\\\^Q 
R^^ w^oo 2R{2R + ^e[-w,w]'^ 



(3) 



proves our lemma. Let us consider first the case s = in the above equation. 

1 



2R{2R 



xe[-w,VF]'* 



1 



2R{2R+l) 
1 



xei-w.VF]'* 



which is obtained by Jacobi identity can be rewritten as, by letting — > oo, 

1 



(4) 



w^cS^^ 2R{2R 



Twri { [Hr+1, -i ^'(x, X + ei)] 

^ I {^,:x.+ei}cl-R-l,R+l]d\l-R,R]'i 



-[Hr+i - Hr, -i ^'(x, X + ei)] } , 

{x,x+ei}c[-R-l,ii+l]'* 

whose norm is estimated as 0{1/R) and goes to zero as i? — >• oo. Next we consider the 
case of finite s in (3). Thanks to 



1 



2R{2R+1Y 
1 



H Xin{x),-ias{[HR,HR+i 

:x.e[-W,W]d 



2R{2R 



xe[-w,w]<* 



< 



+ 



1 



2R{2R+l) 
1 

2R{2R+lY-^ 



' yi£[-W,W]<^ 



(ia„(n(x)) 
du 



, —i[HR, Hr+i] 



to show that the term 
1 



2R{2R 



iir . dau{n{^)) 

^1 xe[-vy,w]<* 



(5) 



for each u converges to zero as — > oo and i? — > oo proves the claim. Since ^""^"^^^^ = 
Ef=i(Q;«(j(x - e;, x)) - Q;„(j(x, x + e/))) holds, it follows that 



(iQ;„(n(x)) 
du 



yL&[-W,WY 
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w-i 

- Wau{3{W,X2,---,Xd,W + l,X2,---,Xd))+ X] x + ei))} 

xi=-W 

d 

+ E E 

1=2 — W<xi,X2,---,xi_i,xi^i,---,Xii<W 

xi {au{j{xi, X2, ■ ■ ■ , xi-i, -W - 1, xi+i, ■ ■ ■ ,Xd,Xi, - ■ ■ , -W, ■ ■ ■ , x^)) 

- au{j{xi,X2,- ■ ■ ,xi-i,W,xi+i,- ■ ■ ,Xd,xi,- ■ ■ ,W + 1,- ■ ■ ,Xd))} (6) 

holds. We employ repeatedly the following group velocity lemma (Theorem 6.2.11 of [11]). 
There exists a positive constant V such that for strictly local observables A and B the 
following inequality holds: 

\\[A,at o T^{B)]\\ < C\\A\\\\B\\ cxp (-|a;| + V\t\) , 

where C is a constant depending only on size of the regions where A and B live in. Let us 
substitute (6) for (5) and estimate its each term. Thanks to || — i[HR, Hji+i]\\ = 0{R'^~^), 

2i?(2i?+ 1)^^-1 II U<.^,.,<,, 

-WauUi-W - 1, xa, ■ ■ • , -W, X2,---, xa)), -i[HR, 



< O (^^0{W'^)ex.-p{-{W - R) + V\u\) 



holds and it goes to zero as 1^ — > oo. The commutator related with 
—Wau{j{W, X2, - ■ ■ ,Xd,W + 1,X2, - ■ • , Xd)) also vanishes in the same manner. Next one 
can see that 



w-i 



2D/2i? + iy-i ll^ ^ ^ <^u{j{^,^ + ei)),-i[HR,HR+i] 

V / X\=—W —W<X2,---,Xd<W 

, (s.-'^')e-o(i 



M 

holds and it goes to zero as R ^ oo. Finally, 

I ||[ V 

2R(2R+ l)d~i^^^ ^ 

V / —W<xi,X2,---:Xi_i,xi^i,---,Xd<W 

Xiau{j{xi,X2, ■ ■ ■ , xi_i, -W - 1, xi+i, ■ ■ ■ ,Xd,Xi, - ■ ■ , -W, Xd)), 

^d\-{W-R)+V\u\ 



holds and it goes to zero as 1^ ^ oo and the term related with xi — W vanishes similarly. 
Thus we proved the equation, 



tL^^oo 2R(2R + ir'^ ^^-' xM^),HR{t)]) = me^)). 

-r ±j xe[-W,W]'' 
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Q.E.D. 

By using an expression Hji — ^^.^i Y^-R<y^<R-i '^-R<yi,-,yr-i,yr+i,-,yd<R 

h{y,y + er), 

one obtains for translationally invariant states, 

-i ^ 

{J{0,e^)) = lim— ^— --^^^ ^ ^ 

z/ti^z/t -h 1) ^^^^ ^^-^ -R<yr<R-l -R<yi,-,yr-i,yr+i,-,yd<R 

xi{[n{x.),at{h{y,y + er))]) 

from lemma 1. Thanks to the translational invariance of the state, substitution z := y — x 
makes the above equation 

-i ^ 

(i(0,ei)) = ^im— ^— — E 

z-ft^z-a i; ^g2;d r-=l -R<yr<R-l -R<yi,-,yr-uyr+i,-,yd<R 

iVi - zi){[n{0),at{h{z,z + er))]). 
One can perform summation for y of the above equation. For r — 1, one obtains 

and for r ^ 1, 

2R{2rI lY-^ ^ ^ (yl-^l)=^^l■ 

zjil^zji Lj -R<y^<R-i -R<yu-,yr-i,yr+i,-,yd<R 

Thus we obtain 

0-(0, ei)) = i ^ I + ^) ([n(0), ai(/i(z, z + ei))]) 

+ ^ E ^i([^(0)> «t(^(z> z + er-))])} ■ (7) 

r=2 J 

Thanks to this representation, we obtain the following theorem. 

Theorem 1 For translationally invariant states with nonvanishing current, there exists 
at least one r e {1, 2, • • • , d} such that neither {[n{0), at{h{z, z + e^))]) nor a dynamical 
correlation function 

{n{0)at{h{z, z + e,)))'^ := {n{0)at{h{z, z + e^))) - (n(0)) (at(/i(z, z + e,))) 

are absolutely integrable with respect to z and t. 

Proof: First let us note that if {[n{0),at{h{z,z + e^))]) is not absolutely integrable, 
{n{0)at{h{z, z + e^)))^ is not also absolutely integrable. In fact thanks to 



{n{0)at{h{z, z + e,)))^ - {n{0)at{h{z, z + e,)))^ = ([n(0), at{h{z, z + e,))]). 
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I ([n(0), Q;t(/i(z, z + e^))])! < 2|(n(0)at(/i(z,z + e^)))^! is obtained. Now we show that 
there exists r such that {[n{0),at{h{z,z + e^))]) is not absolutely integrable. Prom (7), 
one can see that 



l(j(0,ei))| < E{ 



1 



[n{0),at{h{z,z + e,))])\ 



r=2 ) 

< E (ki| + ^)i:i(KO),a,(Mz,z + e,))])| 

holds. Thanks to the group velocity lemma, there exists a constant C and V such that 

|([n(0), at(h(z, z + e,))])| < C||n(0)|| \\h(0, + e,)|| exp (-|z| + V\t\) 
holds. Therefore, for an arbitrary e > 0, there exists N > satisfying 



E 

|z|>y|t|+Ar 



Zi\ 



J2\{[niO),at{h{z,z + 



< e. 



r=l 



Thus we obtain 



IO-(o,eO)| < E (l^il + |)El(KO),«t(Mz,z + e.))])| 

zeZd ^ ^ r=l 

< e+{v\t\+N + l]Y: E |([n(0),a,(Mz,z + e,))])| 

r=l\z\<V\t\+N 



< e+(v\t\+N+l)j2Il |(KO),ai(/i(z,z + e 



which leads to 



i:EI(Ko),«*(Mz,z + e.))])|> ^ 

r=i zezd 



+ iV+ 1/2 



If for all r = 1,2, |([n(0), Q;t(/i(z, z + e^))]) | are integrable, their summation 



Sr=i z + Gr))])! must be integrable. It ends the proof. 



Q.E.D. 



4 Conclusion and Outlook 

In this paper, we considered states over rf-dimensional infinite spin systems which are 
translationally invariant and have non-vanishing current expectations. The dynamical 
correlation function between n{0) and at{h{'x, x + e^)) for some r shows slow decay with 
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respect to space and time (x, t). The key observation to prove our theorem is to ex- 
press the current operator by a commutator and to apply the argument of [14] for the 
continuous symmetry breaking case to our nonequihbrium case, although no continuous 
symmetry breaking takes place in our case. The behavior of dynamical correlation func- 
tion is related with observable effect through response functions, which will be discussed 
elsewhere. Since the conditions for our nonequihbrium states are weak, our result is gen- 
eral but is not very strong. It is interesting to investigate whether additional physical 
conditions can derive more detailed form of the correlation function. 
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